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Abstract 

We consider a description of propagators for particle resonances which takes into ac- 
count the quantum mechanical interference due to the width of two or more nearby 
states that have common decay channels, by incorporating the effects arising from 
the imaginary parts of the one-loop self-energies. Depending on the couplings to the 
common decay channels, the interference effect, not taken into account in the usual 
Breit-Wigner approximation, can significantly modify the cross section or make the 
more long-lived resonance narrower. We give few examples of New Physics models for 
which the effect is sizable, namely a generic two and multiple Higgs model and neutral 
vector resonances in Higgsless models. Based on these results we suggest the imple- 
mentation of a proper treatment of nearby resonances into Monte Carlo generators. 



1 Introduction 



The Breit-Wigner (BW) approach to resonances in particle physics allows to take into ac- 
count the finite width of a meta-stable particle. Various forms of this approach exist, which 
allow to describe different situations, from the narrow width approximation, up to broad 
and energy dependent widths. 

In the following we consider a generalisation of the Breit-Wigner description [1] which 
makes use of a matrix propagator including non-diagonal width terms in order to describe 
physical examples in which these effects are relevant. Indeed for more than one meta-stable 
state coupled to the same particles, loop effects will generate mixings for the masses as 
well as mixed contributions for the widths (imaginary parts). In general a diagonalisation 
procedure for the masses (mass eigenstates) will leave non-diagonal terms for the widths. 

Usually non-diagonal width terms are discarded. This is a good approximation in most 
cases and different unstable particles are described each by an independent Breit-Wigner 
profile. However when two or more resonances are close-by and have common decay channels 
(a precise definition of nearby resonances will be given in the following) such a description 
is not accurate anymore. Indeed, already from a quantum mechanical point of view, one 
cannot treat these states as independent. The usual Breit-Wigner approximation amounts 
to sum the modulus square of the various amplitudes neglecting the interference terms. When 
there are common decay channels and the widths of the unstable particles are of the same 
order of the mass splitting, the interference terms may be non-negligible. Generalisations of 
the Breit-Wigner approach were discussed in the literature in the past (see [1] for a general 
discussion of unstable-particle mixing, gauge invariancc issues, unitarity and renormalisation 
for strongly-mixed systems), mainly focussing to the applications in the Higgs sector of the 
supersymmetric standard model and CP violation [2-6]. In the following we shall give 
a general formalism for scalar and vector resonances and discuss other relevant physical 
examples. In particular we will consider models of physics Beyond the Standard Model 
(BSM) in which new resonances play a crucial role, such as a generic two-Higgs model or 
the case of an arbitrary number of Higgs particles. Interesting is also the case of Higgsless 
models for which the lowest lying Kaluza-Klein (KK) excitations of the photon and the Z 
are nearly degenerate and have common decay channels into fermions and gauge bosons. We 
will show that in all these cases, the interference effects can play an important role. Based 
on these results we suggest that a proper treatment should be carefully implemented into 
Monte Carlo generators as physical results may be dramatically different from a naive use 
of the Breit-Wigner approximation. 

2 Scalar fields 

We first discuss the case of scalar fields which gives a simpler overview of the problem 
without the extra complications of the gauge and Lorentz structure present in other cases. 
The main observation needed to take the mixing effect into account, is to remember the link 
between the action and the propagator in quantum field theory and write the correction to 
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the propagator as a modification of the kinetic operator, which, in the general case, is a 
matrix. 

Let's start with the action for a real scalar field given by (in position and momentum 
space) : 



-'scalar 



I d'x \ [(9,0)^ - ml<^^\ = / \4>{-v)KMm , (2.1) 



where, for later convenience, we have defined a kinetic function 

K^M^p'-ml (2.2) 
The propagator of the scalar field is defined as the inverse of the kinetic operator: 



iA,,o = iK-'o = ■ (2-3) 

In the presence of interactions, the kinetic term will receive contributions by loops: if we call 
in(p^) the value of the 1-PI corrections to the propagator, the corrected kinetic term is 

i^,=p2-m^ + n(p2). (2.4) 

Therefore the new propagator is 

^A. = ^K;^ = .2 _ = ^^.o ("A^o)'' , (2.5) 



p2 — 777.0 + n(p2) 



n=0 



which corresponds to the resummation of the 1-PI insertions on the bare propagator. The 
pole of the resummed propagator defines, as usual, the renormalised mass. If the particle is 
unstable, n(p^) is complex. The real part is used to renormalise the mass and the imaginary 
part defines the width of the particle. The propagator has a pole in — imT, with m the 
renormalised mass and F the width. For our purposes, we leave out the imaginary part 

QU(p') = E(p^) . (2.6) 

In the narrow width approximation, E(m^) = mF. 

This formalism can be generalised to a system involving multi-fields, which do couple 
to the same intermediate particles: the loops will generate mixings in the masses, but also 
out- of- diagonal imaginary parts. In general the real and imaginary parts will not be diago- 
nalisable at the same time: we are interested in the phenomenological consequences of this 
scenario, when the out-of-diagonal imaginary part is of the same order as the mass sphtting. 
The kinetic function is now written in matrix form: 

{Ks)ik = (p' - m^)6ik + ^^lk{p') ■ (2.7) 

(We are considering the imaginary part only, the real one is used to renormalise the masses.) 
The propagator of the fields can be defined as the inverse of this matrix: 

i(As)ik = i {k;')^^ . (2.8) 
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To give an explicit example, let us focus on the two-particle case: 



22 

i^2i P 



(2.9) 



where 



Ds = (p' -mj + iEn)(p' -ml + iEss) + S12E21 . 



(2.10) 



For vanishing S12 and S21, the propagator is diagonal and it reduces to two independent 
Breit-Wigner propagators with miTj = Sjj(m^). 

However, the narrow width approximation is not valid if the off-diagonal terms are sizable 
compared with the mass splitting. Defining 2M^ — ml + ml and 25 = m| — mf , the poles 
of the propagator (zeros of Ds), which define the physical masses and widths of the two 
resonances, are: 



m 



- i5ii±^ ± ^V(S22-Sn + 2^5)2 + 4E12S21 



(2.11) 



Note that the value of the masses is modified by the presence of the off-diagonal terms due 
to the imaginary part of the square root, at the same time the widths are affected. More 
importantly, the off-diagonal terms in the propagator will generate non-negligible interfer- 
ence, which can be in turn constructive or destructive. This effect will be illustrated with 
some numerical examples. 

Finally, let us write some general formulae for S(p^) . If we assume that the scalar particles 
i and j couple to a pair of particles a with couplings and A;^ respectively, the matrix Ej^ 
can be written in general as 



(2.12) 



We will consider couplings with fermions /, scalars s and vectors V: 



(2.13) 



where Pl,r = (1 T 75)/2 are the chirality projectors. Therefore, expanding for small masses 
of the particles in the loop (full results are given in appendix A), we get: 



e(;)(p^) 



167r 



+ ... 



IGtt 
647r rriy^my^ 



+ ... 



(2.14) 
(2.15) 
(2.16) 
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Note that the momentum dependence in the previous formula could give rise to violation 
of high-energy unitarity at energies above the resonance as well as distortion of the line-shape 
for broad resonances. These issues are discussed in detail in [1,7]. Actually we do not face 
these problems here, since we are interested in the effects of the absorptive self-energies near 
the pole of nearby resonances whose widths are of the same order of magnitude of the mass 
splitting. 



2.1 A numerical example: near-degenerate Higgses 

As a numerical example, we will study two heavy Higgses where both the scalars develop a 
vacuum expectation value (VEV) and therefore couple to the W and Z gauge bosons. This 
situation is common in supersymmetric models where two Higgses are required by writing 
supersymmetric Yukawa interactions for up and down type fermions, and generic two Higgs 
models. The interference between near degenerate Higgses has been studied in [2-4] focusing 
in CP violation effects. 

The couplings of the two CP-even Higgses to gauge bosons can be written as 

^wwHi = 9 rnw cos a , )^wwh2 = 9 rnw sin a , 

AzzHi = —cos a, \zzH2 = sma; (2.17) 

cos Uw cos Uw 

where a is a mixing angle taking into account the mixing between the two mass eigenstates 
and the difference between the two VEVs. 

Here we are interested in a generic production cross section of the two nearby Higgses 
on the resonances, with decay of the Higgses into gauge bosons (either WW or ZZ). The 
amplitude of this process is proportional to the resonant propagator weighted by the cou- 
plings given in eq. fl2.17p . In the case we are considering, the common decay channels can give 
off-diagonal terms in eq. (12.91) which are sizable compared with the mass splitting. Therefore 
we need to include their effects and a generic cross section will be proportional to (here we 
assume that the coupling to the initial particles are the same): 

|(A^i + ^f) cosa + (Af + sina|' . (2.18) 

In Fig. [H we plot this quantity in arbitrary units and compare it with the Breit-Wigner 
approximation: we fix mni = 400 GeV, and vary the splitting from 50 to 5 GeV. For 
simplicity, in the following we will assume a = 7r/4, so that the two scalars have the same 
couplings (but this assumption is not crucial for our conclusions) . The exact treatment of the 
resonances unveils a destructive interference (which is neglected in the naive Breit-Wigner 
case) that can drastically reduce the cross section. Also, the interference between the two 
resonances splits the mass poles [3,4]. 

This effect can be even more important for scenarios with a large number of scalars 
as predicted in some string models. Our analysis can be easily extended to an arbitrary 
number of Higgses. Let's take for example the couplings to the gauge bosons to be given by 
5'sm/ V^, where 9sm is the SM coupling of the gauge bosons and is the number of Higgses. 
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Figure 1: Plots of the production cross section (in arbitrary units) of two nearby Higgses decaying into 
gauge boson pairs for the naive Breit-Wigner (blue-dashed) and exact mixing (red-solid). The mass of the 
first resonance is fixed to 400 GeV, the sphtting respectively 50, 25, 10 and 5 GeV and a = 7r/4. 



In Fig. [2] we plot the cross section for seven nearby Higgses, with the first one at 400 GeV 
and the others at a distance of 5 and 10 GeV, the width of each being 6.2 GeV. From the 
plot it is clear that the destructive interference reduces the giant resonance (which is not 
distinguishable from a single Higgs, once the experimental smearing is taken into account) to 
a bunch of gnometti (dwarfs), which will be very hard to detect. The cross section is in fact 
reduced by a significant factor with respect to the naive expectation, and the smearing will 
wash out the peak structure. This situation is different from the continuum Higgs spectrum 
of [8], where the Higgs peaks are smeared by the experimental uncertainties only. Therefore, 
in this case, together with the appearance of a "continuum" , the cross section is suppressed 
by the interference. It is intriguing to compare this analysis with Un-Higgs models [9,10], 
where the Higgs in indeed a continuum: such behaviour may arise from the superposition of 
Kaluza-Klein resonances in extra dimensional realisations or deconstructed models [10]. 



3 Vector fields 

For a vector field the technique is similar to the scalar case, however a major issue is defining 
the propagator of a metastable particle in a gauge invariant way (see for example [11]). 
Nevertheless, at the pole, only the values of the pole mass and width (which are gauge- 
invariant) are relevant, and extra terms that must be added to preserve gauge invariance are 
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Figure 2: Plots of the production cross section (in arbitrary units) for seven nearby Higgses equally coupled 
to SM gauge bosons: the naive Breit-Wigner (blue-dashed) bump reduces to a row of seven dwarfs when 
the exact mixing (red-solid) is taken into account. The mass of the first resonance is fixed to 400 GeV, the 
sphtting between the six Higgses respectively 10 and 5 GeV. 



numerically negligible. The kinetic function for a vector V, in generic ^-gauge, is 
K^^^o = Qi^uip^ - ^o) - Pf^Pu (7-1 



^.-^)(/-Mo^) + ^(f-Mo^)- (3.1) 
The propagator, is then defined as the inverse of the kinetic term: 

Pf,Pu\ -I Pi^Pu -i^ . . 

The first term of the propagator has a gauge-independent pole at the V mass, while the 
other part has a gauge-dependent pole, which will cancel the pole given by the Goldstone 
boson (whose mass is indeed ^Mq). Therefore, at the pole we can neglect the contribution 
of the second term, and the propagator simplifies to 

p2 J p2_ Ml 



iA^,,o -[g^^u-^] o _ . .2 ; (3-3) 



which is equal to a Lorentz tensor times a scalar propagator. 

Factorising out the Lorentz structure, loop corrections can be parameterised as 

= Ut gf,u + p^p^ , (3.4) 

so that 

K,.,o + n,. = (^g,, - ^) ip' - Ml + + Ml + + p^n^^ . (3.5) 



The corrected propagator is therefore: 



iA 



+ 



(3.6) 



The first term defines the pole mass and width (gauge independent), while the second term 
contains a gauge-dependent pole, and it is negligible at the physical pole. Here we will be 
interested only in the imaginary part = which defines the decay width of the vectors. 
Neglecting the second part of the propagator, we find: 



iA 



P^ 

p2 



(-2)A,(My). (3.7) 



For a generic number of vectors, the same discussion as in the scalar case applies, and 
has a matrix form, like in (12. 8p . 

We conclude with some general formulae for S which can be expressed as in eq. (I2.12p . 
We will consider couplings with fermions /, scalars s and vectors V: 



fiup / ; 



(3.8) 



where Gf,„p = g^y{p + qi)p + gup{q2 - qi)fi - QpM^ + P)u, with p = q^ + ga- Expanding for 
small masses of the particles in the loop (full results are given in Appendix A) we obtain: 



J I] 



247r ^ 



[^'PiP%, = ^P' + --- 



487r- 

^v^v P 
1927r mlrm\ 



+ 



(3.9) 
(3.10) 
(3.11) 



In the following we apply this formalism to a simple case. 



3.1 Numerical example: Z' and A' in Higgsless models 

In Higgsless models [12,13], the first two neutral resonances are nearly degenerate, and they 
correspond to the first KK excitation of the Z and of the photon. Here we will explicitly 
refer to the warped extra dimensional model in [13]: the masses can be approximated by 

m|, ~ m^;^ + 4m| , m^, ~ m^-j^ , (3.12) 

so that the mass difference is very small: 

/lTeV\^ 

mz'-mA'^2 — ^~16GeV- . (3.13) 

rriKK \ rriKK J 
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In terms of the parameters of the warped geometry {R is the curvature, R' the position of 
the Infra- Red brane in covariant coordinates): 

2 4 1 

mxK ~ ^ , mw = — jfT ; (3.14) 

R R' log 

therefore, given the value of the curvature R, the KK mass {R') is determined by the W 
mass. 

The coupling to the light gauge bosons can be estimated by use of the Higgsless sum 
rules, that ensure the cancellation of the terms growing like the energy square in the elastic 
scattering amplitude of the longitudinal W and Z, thus delaying the violation of perturbative 
unitarity at higher scales than in the Standard Model without Higgs. One of these sum rules 
is [14,15]: 

gwwww = -7 i 9wwz—i- + 9wwk—r ] + ^ ' {3.15) 

where gwwww = 9^, 9wwz = g cos 9w and we have included a partial contribution for the 
Higgs (C = 1 corresponds to the Higgsless limit, C — ^ to the SM Higgs). Neglecting the 
contribution of the heavier states, from eg. (13. 151) we can estimate the coupling of the first 
tier: 

9wwi-\h:9 • (3.16) 

V 3 niKK 

This is actually the coupling of the first KK mode of the neutral component in the SU(2) 
multiplet: therefore, in order to obtain the couplings of the mass eigenstates, one needs to 
impose a rotation by an angle 6'i, which describes the mixing between the SU(2) and the 
U(l) vectorial component: 

gwwz' - 9wwi cos 01 , QwwA' - 9wwi sin 6i . (3.17) 

Numerically, it turns out that the Z' is 45% in the SU(2) and 55% in U(l), therefore cos^i ~ 
v/45% ~ 0.67 and sin^i ~ ^55% ~ 0.74. 

The amplitude for the decay Z', A' — > WW is given by {aw = g'^/in): 

2x f cos^^*! cos^i sin 6^1 \ (aw ^oion 
)-[ cos^i sin 01 sin^^^i J 144 ml^mlj, ' ^ ^ 

and the width can be estimated by: 

^ S(m2,^) f 17 19 \ ^ fmKK\^ 



rriKK 



19 21 ■ {rfeVJ 



GeV- C- (3.19) 



The off-diagonal entries are large, and they are of the same order of the mass splitting 
between the two neutral gauge bosons. 
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R 


mw' tua' tuz' 


gwWA' gwWZ' QffA' 9ffZ' 


10-^ 


1032 1028.5 1043.79 
805.3 801.0 820.9 
634.6 629.0 655.0 


0.022 0.019 -0.021 0.065 
0.028 0.024 -0.027 0.085 
0.037 0.029 -0.035 0.11 



Table 1: Three Higgsless points for the model proposed in [12], masses in GeV. 



The determination of the couplings to fermions, which are relevant for the Drell-Yan 
production, is more involved. In fact, the typical scenario is that the right-handed com- 
ponents of the light quarks and leptons are locahsed on the Ultra- Violet brane, while the 
left-handed components are spread in the bulk [13]. For the left-handed couplings, we have 
some freedom: however, electroweak precision measurements prefer small couplings therefore 
we will neglect them. The couplings of the right-handed components depend uniquely on 
the suppression of the wave functions on the UV brane. Therefore, we can approximate: 



9f!i -Qf9 tanOw 



Numerically we find that: 



QfQ tanOw 2.4^^ ~ Qf ■ 0.066 ■ ( . 



9ffZ' - 9ffi 



rriKK 



9ff-A' = 9ff-iV ; 



(3.20) 



(3.21) 



were the ratio r) ~ —0.32 is to a good approximation independent on the value of the KK 
mass. In Table 1, three numerical examples are given, corresponding to different values of 
the curvature: the masses and couplings are calculated exactly, following [12], and confirm 
our estimates. 

We are interested in three processes: Drell-Yan production and decay into gauge bosons 
W'^W~ (DY), Drell-Yan production and decay into a pair of leptons (Leptonic) and vector 
boson fusion production followed by decay into gauge bosons (VBF). As in the scalar case, 
the amplitudes of these processes on resonance are proportional to the propagators weighted 
by the couplings with the incoming and outcoming particles. The cross section a{qq — > 
A', Z' — > Vt^+Vt^") is proportional to: 



cos6li Af ^' -Fsin6liA 



Z'A' 
s 



?7|(sin^iAf^'+cos^i Af^') 



A'Z'^ 



which must be compared to the naive sum of two Breit-Wigner distributions: 

2 



'DY 



cos 9i A 



Z'Z' 



\r]\ sin 9i A 



A'A' 



(3.22) 



(3.23) 



On the other hand, the leptonic cross section a{qq — > A', Z' —>■ ) is proportional to: 

2 



CLepto 



Leptonic 



r^^Af- +A 
T)' Af ^' + A, 



Z'Z' 

s 

Z'Z' 
s 



ri\ (Af^'+Af-) 



A!Z' 



(3.24) 
(3.25) 
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Figure 3: Plots of production cross section (in arbitrary units) of the two low- lying neutral resonances 
of the Higgsless model for the naive Breit-Wigner (blue-dashed) and exact mixing (red-solid). The rows 
correspond (from top to bottom) to DY, Leptonic and VBF; the columns (from left to right) correspond to 
mKK = 1000 GeV, 800 GeV and 600 GeV. 



Finally, for the VBF channel we have: 



O'VBF 



naive 
^VBF 



sin'e, Af^ +cos^^^i Af ^ +cos^^i sin^^i (Af ^ + A 

2 



A'Z' 

s 



sin' 01 Af^'+cos^^i A 



Z'Z' 

s 



(3.26) 
(3.27) 



DY and Leptonic cross sections are easily calculable; the VBF channel requires a numerical 
study or a more involved approximate analytical expression. 

In Figure [3] we plot, for illustrative purposes, the squared matrix element of the three 
resonant production channels for A' and Z' as function of ^/s for three different cases: 
f^KK = 1000 GeV, 800 GeV and 600 GeV. For large masses, the effect of the interference 
is very important and it can affect the value of the cross section significantly. To make this 
more clear, we give hereafter the ratio of the area under the peaks in the figure obtained by 
the exact formula and the BW case. This roughly corresponds to the ratio of the integrated 
cross sections. 
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Mkk^ 1000 GeV 800 GeV 600 GeV 



DY: 



1.6 

3.15 
0.6 



1.15 
1.4 
0.8 



1.02 
1.05 
0.97 



Lept: 



VBF: 



In the VBF channel there can be a reduction up to 50%, while in the other two channels 
the interference is constructive and the total cross section can be enhanced by a factor of 
2-3. The interference is therefore extremely important, especially in the TeV region. Since 
this represents the upper bound for Higgsless models, the interference effects are crucial to 
determine if the whole Higgsless parameter space can be probed at the LHC. As a conse- 
quence the implementation of the exact matrix propagator in the Monte Carlo generators 
seems mandatory for a correct analysis. 



We have shown that for two or more unstable particles, when there are common decay chan- 
nels and the masses are nearby, the interference terms may be non-negligible. This fact is 
well known in the literature and well studied since more than a decade [1], here we provided 
further examples in which this formalism is important. This kind of scenario is not uncom- 
mon in models of New Physics beyond the Standard Model, especially in models of dynamical 
electroweak symmetry breaking or in extended Higgs sectors. We reviewed the formalism 
for scalar and vector fields based on a matricial form of the propagator for multi-particles 
taking care of this case. This formalism can be easily extended to fermionic resonances as 
discussed in [1] and applied to heavy neutrino mixing in [16]. We gave few examples in which 
the effect of the non-diagonal width is important in physical results. In models with multi- 
Higgses and in Higgsless models with near degenerate neutral vector resonances, we showed 
that interference induced by the off-diagonal propagators are very important and they can 
either suppress or enhance the total cross sections on resonance depending on the relative 
sign of the couplings to the initial and final states. Similar issues were discussed in the CP 
violating MSSM for the Higgs sector in [17] and [18]. We expect similar effects in all the 
New Physics schemes involving resonances for which the mass sphtting is of the same order 
of the off-diagonal matrix elements in the propagator due to common decay channels. For 
example in generic Technicolour models in which vector and axial-vector composite states 
can have the aforementioned property or also in supersymmetric models with two nearby 
neutralinos. Other examples are models in warped extra dimension, like gauge-phobic Higgs 
models and Composite Higgs models. 

As a conclusion, the interference effects can be crucial to study the phenomenology of 
such models at the LHC, and to determine its discovery potential. We stress that a proper 
treatment should be carefully and systematically implemented into Monte Carlo generators 
used to study BSM models, as physical results may be dramatically different from a naive 
use of the Breit-Wigner approximation. 



4 Conclusions 
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A Appendix: exact amplitudes 

We give in this appendix some detailed formulae which were not given or given in approximate 
form in the text. Defining: 



X{p,mA,mB) = 



(A.l) 



the imaginary contribution to of the 1-PI corrections of the scalar amplitudes are, at one-loop 
level, {mA and ms are the masses of the particles in the loop): 



■s/X{p, mA,mB) 



IGtt 

-2iX^Xi + X^^)mArnB] , 

^/X{p, mA^m s) 
IGtt 

yJX{p, mA,mB) 



[(AlAi, + A^«Ai)(p^-m^-m^)+ 



-XlXi 



647r 



Ay Ay 



m 



2 \2 
B) 



+ 



In a similar way for the vector amplitudes we have: 



i/A(p,mi,m2) 



AlAI + A^Ajj / _2 _2 



IGtt 

+ (AiAij + A^Ai)2m^mB] , 

^A(p, mA^ms) X\X{ 
IGtt 3 

{X{p,mA,mB)f''^ 



2p — — 



6 



rriA — m 



2 \2 



P" 



1927r 



^V^V 2 2 



m\ + m% m\ + m% + I0m\m% 



p-* 



(A.2) 
(A.3) 

(A.4) 



+ 



(A.5) 
(A.G) 



(A.7) 
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